Hubbard model on d-dimensional hypercubes: Exact solution for the two-electron case 
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The Hubbard model is exactly solved for two particles with opposite spins on d-dimensional 
hypercubes. It is shown that the spectrum can be separated into two parts: a trivial (U-independent) 
part resulting from symmetries of hypercubes and a non-trivial part expressed as a single-impurity 
problem on a set of finite chains of size d' + 1 {d' < d) . The exact expression for the one-particle 
Green's function is given. Finally, we discuss the extension of these results to standard hypercubic 
lattices with periodic boundary conditions. 
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One of the most widely used models to describe 
strongly correlated fermions systems is the Hubbard 
model and its various extensions. Unfortunately, only 
limited exact information about its physical properties is 

' available. For the one-dimensional lattice the celebrated 
solution of Lieb and Wu ^ provides the exact eigenspec- 
trum of the model. In the opposite limit of large dimen- 
sions investigated very recently some almost exact 
results have also been obtained. However, for intermedi- 
ate dimensions, and particularly for the very important 
two-dimensional case believed to be relevant to high-Tc 

] superconductivity, very little is at our disposal. Most 
of the results reported so far have been obtained either 
from numerical solutions on very small clusters (subject 
to important finite-size effects) or by using a variety of 
approximate analytical methods (with domains of valid- 
ity and/or systematic errors difficult to evaluate) 

In this paper we present an exact solution for the Hub- 
bard model with first-neighbor hopping term t and on- 
site interaction energy U 

: ^ = E + E 4^*t4c'U = T{t) + v{u), 

, <ij><y i 

: (1) 

for two particles with opposite spins on a d-dimcnsional 
hypercube 7^. The hypercube 7^ is defined as the set 
of iV = 2'^ sites whose d coordinates are either or 1. 

\ Although our solution is yet limited to the two-electron 
case, to exhibit exact results for a truly interacting sys- 

[ tern in dimension greater than one is clearly of primary 
interest. This is particularly true since hypercubes are re- 
lated to usual cubic lattices with periodic boundary con- 
ditions: 72d is topologically equivalent to a d-dimcnsional 
hypercubic lattice of linear size equal to 4 with periodic 
boundary conditions, noted Zf (i.e., 72 is equivalent to 
the 4-sites ring, 74 equivalent to the 2-dimensional 4x4 
cluster with periodic conditions 0,^, etc.), while ^2d+i is 
equivalent to two translated linked copies of Zf. To solve 
the above Hamiltonian we first use an abelian subgroup 



of the full 7d point group, to block-diagonalize the initial 
matrix of size 2^"* into 2'^ blocks of size 2'^. Quite interest- 
ingly, these smaller submatrices correspond to a simple 
family of effective one-electron Hamiltonians defined on 
Id' {d! < d), some of which having a single-impurity site. 
In a second step, we show how to further reduce these 
matrices into smaller blocks of size d' -I- 1 (instead of 
2''), corresponding to finite chains of size d' with one 
impurity site at one end and new specific hopping terms. 
Using a standard approach for impurity problems we pro- 
vide a closed expression sum rule for the eigenspectrum. 
We also derive the exact expression for the one-particle 
Green's function. Finally, the extension of these results 
to hypercubic lattices with periodic boundary conditions 
is briefly discussed. 

Let us define A^;, an abelian subgroup of the full jd 
point group, generated by the d reflections into perpen- 
dicular (d — 1)— dimensional planes meeting at the 7^ 
centre. 

Ad = {TTi, i = 1 . . . d I TTi^ = TTiTTj = TTjTTi } (2) 

where E is the identity operator. It is easy to show that 
Ad has 2'^ elements (owing to the commutativity between 
the mirror operations), and that the orbit of a generic 
point under A^ is a hypercube 7^. Let us number the 
Ad elements L by integers between and 2^* — 1, in the 
following way: L being a product of mirrors tt^, let L 
be the number, written in base-two, whose correspond- 
ing i*'^ digits are equal to one (tti is L = 1, 7ri7r27r3 is 
L — 7, etc.). Let us also locate the sites according to 
their numbering in base-two: The i*'* coordinate is equal 
to the i*^ digit. Now, the action of the symmetry op- 
eration L, onto a site s, noted as the "special product" 
L • s, can be explicited: each time a digit equals 1 in 
L, it switches the corresponding digit in s. For example, 
in 3 dimensions, 7 • 1 = 6, with corresponds in coordi- 
nates to 7ri7r27r3(0, 0, 1) = (1, 1, 0). To construct the irre- 
ducible representations (irreps) and the character table 
of Ad, it is useful to remark that this group is isomor- 
phic to d times the tensorial product with itself of the 



1 



two-element group Zi = {i5,7ri} = Ai ||^. This is noth- 
ing else than the group theory counterpart of the natural 
construction of 7^ as two displaced copies of 7d-i- The 
Acj irreps are easily obtained as tensor products of the 
two irreps Fq and Fi of Zi. As an illustration, the char- 
acter table fo A2 = ^2 X Zi is shown in Table |. All 
the irreps are uni-dimensional (A^ is abelian), and the 
character table is identical to standard Hadamard matri- 
ces 0, which are orthogonal matrices with elements ±1 
(note also that lines and columns of Hadamard matrices 
represent the non-interacting one-electron orbitals of the 
problem). Note that the iV (= 2'*) irreps can also be 
labelled by integers M, through their base-two decompo- 
sition, by ordering the different occurences of Zi in the 
tensorial product, and fixing the corresponding digit to 
or 1, whether the irreps Fq and Fi occur. 



TABLE I. Character table for A2. The irreducible repre- 
sentations arc numbered by an integer M as defined in text. 
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Let us consider first the one-electron tight-binding 
spectrum with hopping term t. Each eigenstate belongs 
to a given irrep of A^, and reads 

N-l 

l^'^>=7P E x^^i?(|o» = -=E(-in''1-o> 

(3) 

where mM,r is the number of digits equal to 1 simul- 
taneously in M and r (we have used the above defined 
special product • to denote the action of the symme- 
try operation R onto the basis kets). Since all the sites 
are equivalent, it is sufficient to consider the interaction 
of the site number with its d first neighbors which 
have exactly one digit equal to 1 and belong to the set 
Vd — {vj = 2^,j = 0. . .d — 1}. As a result, one finds 
that the eigenenergies read 

£^M = 5]x^f = -t(rf-2iVM) (4) 
i=o 

where Nm is the number of digits equal to 1 in M. Since 
M runs from to 2'' — 1, the 7^ spectrum consists of 
d + 1 levels, from —dt to dt, with equal spacing 2t. The 
degeneracy of a Em level is the standard binomial coef- 
ficient counting the number of sets of Nm digits among 
d. As d tends to infinity, one recovers a gaussian-shaped 
spectrum, as for the infinite-dimensional cubic lattice. 

Let us now consider the two-electron case (electrons 
with opposite spins). The full Hamiltonian matrix, of 
size 2^"^, is readily block-diagonalized into 2^^ blocks of 
size 2'', associated with each representation M. For a 
given M, one constructs a basis with 2'^ kets \M, I > 

l^''>=4z E x^'^(|OT,U>) 

^2 ReAd 

1 

= ^J2Xr'\ir'0)^,{r.l)i> (5) 

It is clear that U only occurs in < M,0\V{U)\M,0 >, 
while the kinetic part reads 

< M,l\T\M,l' >= 
^ E E ^r'x^f < (r.O) T (r.O i |T|(/.0) T, (r'.O i> 

r=0 r'=0 
N-l 

= ^ E [X*'xf < (r-O) T (r./) i \T\{r.O) T, (r-/') i> 

d-1 
3=0 
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< (r . 0) T (r • i \T\{vj • r . 0) T, {v, • r • I') l>] (6) 

In the latter expression, we have decoupled the terms 
where the | spin jumps (first part) from those where the 
t spin jumps. In the latter part (sum over j) the only 
nonvanishing term corresponds to r • I — Vj • r • I' —* 
I' = Vj»l, which means / and I' neighbors (recall that the 
group is abelian and its elements are their own inverse). 
One finally gets that the nonvanishing kinetic terms read 



< M,l\T\M,Vj • I >-- 



Ar A.Vj9 



.) (7) 



r=0 



The block matrix associated to the irrep M cor- 
responds therefore to an effective one-electron tight- 
binding Hamiltonian, with sites labelled by I, kets \M, I > 
and hopping terms given by the above expression (^. A 
simple case is provided by the identity irrep M = 0, for 
which = 1 for any r. The effective Hamiltonian iJ'^*'^ 
corresponds to a hypercube 7^, with one "impurity" site 
(with diagonal term U) and a constant hopping term —2t. 
We now show that the other irreps correspond to effective 
Hamiltonian H^^^^ on hypercubes 7^' with d! < d. One 
has to evaluate the right-hand part of equation (^. By 
definition of the Vj (mirrors which connect a site to one of 
its first neighbors), r and Vj • r differ by one digit, which 
depends on j, not on r. If this digit takes the value 1 in 
the base-two decomposition of M, then X^X^j»r = ^1 
for any r, and the corresponding matrix element 

vanishes. H'^^^^ therefore corresponds to a 7^ in which 
families of parallel edges have been cut (those orthogonal 
to the mirrors labelled by the digits 1 in the decompo- 
sition of M). Figure |l| displays graphically this general 
result for d = 2 and d = 3. 
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Fig. 1. Graphical illustration of the effective Hamiltonian 
H^'^^'l associated with the irrep M . Full circles correspond to 
impurity sites with a diagonal term U , full edges to hopping 
terms equal to —2t, and dotted lines to the edges which are 
cut in the given irrep. 

a) d = 2, the four irreps, corresponding to one square with one 
impurity, two times two segments, one containing an impurity 
site, and finally four isolated sites (one being an impurity with 
a diagonal term U). 

b) the case d = 3. Four irreps among the eight are shown. 



So, at this step, the spectrum is solved in terms of 
one-electron tight-binding models on hypercubes of di- 
mension lower or equal to d. The only non-trivial part 
corresponds to those part, in the hypercube decomposi- 
tion, which contain the impurity site We now show 
how to further greatly reduce their complexity. 

To do that it is convenient to re-express the effective 
Hamiltonian H (hereafter we will drop the superscript (0) 
since only the M = case has finally to be considered) 
in the one-electron basis, Eq.(||). It is easy to verify that 
H can be written in the form 



H = Diag{e) + ^In 



(8) 



where Diag{e) is a diagonal matrix whose entries are the 
2'^ free electron-solutions of the hypercube with hopping 
2t (solution at [/ = 0) and In is the N x N matrix with 
unit entries for all i and j. The latter part is nothing 
but a projection operator on the eigenstate \vo > with 
energy U and H can be rewritten 



H = Diag{e) + U\vo >< Vo\ 



(9) 



with {vo)^ = ^ l^l,...,N. 

To proceed further we define for each subspace corre- 
sponding to a given value {d + 1 subspaces of degen- 
eracy gi, the binomial coefficient) a basis consisting of 

the normalized vectors w^"* = —^(1, . . . , 1) and a set of 



vectors {w^*"* , fc 7^ 1} spanning the subspace orthogonal to 

v^^K In this new basis only the (d + 1) vectors vi have a 
non-zero overlap with \vo >, and H, of dimension 2^^, de- 
composes into a diagonal matrix having 2'^—{d+l) trivial 
solutions ei, and a residual U-dependent part, noted 7i, 
given by the matrix of linear size (d + 1) written in the 
form 



n = Diagie) + U\v >< v\ = T^o + V 



(10) 



where \v > has components Vi = ^ gi/N and the e 's 
represent now the (d+1) distinct free-electron energies: 



e^ = -2t{d-2i) i = 0, 



(11) 



An alternative representation consists in going back to 
the basis where the potential operator V is diagonal. H 
is then found to be tridiagonal (one-dimensional tight- 
binding model) with off-diagonal hopping term given by 
ti^i+i = 2t^J [d — -I- 1), i = to d — 1 and a di- 
agonal ^/-contribution at the initial site (here numbered 
0). In other words, the problem is mapped onto a single- 
impurity problem on a finite chain of linear size (d~\- 1). 
Figure |2| illustrates this result for d < 3. 
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Fig. 2. Graphical illustration of the single-impurity one- 
electron problems associated with the non-trivial part of the 
spectrum (general formula given in text). 



It is known that eigenspectrum of single-impurity 
problems can be expressed under the form of a sum rule 
expression using the Koster-Slater approach This 
result can be briefly re-obtained as follows. First, the 
following operator identity is invoked 

u 1 , ,1 



\v >< v\ 



{-Ho - z) 



(12) 



where z is an arbitrary complex number and Gm the 
unperturbed Green's function given by G^^v {z) —< 
v\{T-Lq — z) \v >. Then, using this identity the fully 
interacting Green's function,G^t,(2:), can be written as 



Gvv {z) 



GvJ (z) 



l + UG^"J{z) 



(13) 



Searching for poles of Gm{z) we get 
d 



1 



1 



1=0 



Eq.([l4[) is the final closed expression determining the 
non-trivial part of the spectrum. It is easily shown that 
this equation admits (c? -I- 1) distinct solutions we shall 
denote as SpeS'^\U). 

To summarize, the set of 2^"^ eigenvalues of the hy- 
percube consists of a U-dependent part given by the 
collection of the non-trivial spectra at dimensions lower 
or equal to d: 



e 

i=0 



(15) 



(note that here the binomial coefhcient counts the num- 
ber of occurences of 7^ hypercubes with an impurity site 
in the decomposition of the initial 7^, see Fig.|l|) and a 
trivial U-independent part given by: 



El = 2tl 



I = -d...d 



with degeneracy: 

id 



91 



~td—i 



I 

l+i 
2 



- 1 



(16) 



(17) 



the prime in indicating that the summation over i 
is done with an increment of two. Formula (17) can be 
obtained by tracing back the contributions due to the dif- 
ferent irreps and those issued from the internal symmetry 
of hypercubes (trivial solutions extracted when passing 



from Eq.(^) to Eq.(p^). The ground-state energy is 
not degenerate and corresponds to the lowest solution in 
Spect'^'^\U). Note also that in the limit of large dimen- 
sions, and after proper renormalization of the parameters 
of the Hamiltonian, an infinite-dimensional model can 
be defined, as done recently for a number of correlated 
fermions models, see ||]. 

Let us now consider dynamical properties. We are in- 
terested in evaluating the one-particle Green's function 
defined by 



Gi^\z,U) =<^o\ak 



H 



(18) 



where j-^o > denotes the one-particle ground-state con- 
sisting of one electron of given spin and creates one 
electron of opposite spin in a one-particle state, noted k. 
Hereafter, k varies from to d and labels one of the (d+1) 
degenerate subspaces of the one-electron problem, the 
ground-state corresponding to /c = 0. Now, remark that 
vector a^\tl'o > belongs to the subspace corresponding to 
the decomposition on hypercubes of dimension {d — k). 
There exist 2^ different families of hypercubes having 
such a dimension (See, fig. |l|) among which only one has 
the U-impurity site. Accordingly, we get the following 
result 



(14) Gf{z^U) 



^G 



(d-k) 



{z.U) + {l-^)G, 



[d-k) 



iz,U = 0). 
(19) 



Now, we need to evaluate the fundamental quantity 
G^f \z,U) with Hamiltonian (p^, the dimension d' rang- 
ing from to d. For that we note that the ket |0 >= 
flg li/^o > represents the (first) basis element correspond- 
ing to the diagonal energy eo in representation ( pO| ) . Pro- 
jecting out identity (|l^) onto vector |0 > and expressing 
the different quantities in terms of the non-interacting 
spectrum we get: 



G^"^(z,[/) = 
U 



id') 



ZHI 



2d' 



Ed' 
i=0 ' 



:) 



(20) 



This equation together with Eqs.(pT[) and ( pj| ) give the 
exact one-particle Green's function of the problem. 

Generalization of this approach to higher fillings is 
presently under investigation. We have already found 
that the existence of large fractions of U-independent 
eigenvalues is still valid for some specific fillings. How- 
ever, the underlying structure for the U-dependent part 
of the spectrum is more difficult to elucidate. 

Back to the two-electron case, we would like to men- 
tion that, by using rotations instead of reflections, a sim- 
ilar calculation can be done for standard d-dimcnsional 
hypercubic lattices of linear size L and periodic bound- 
ary conditions {Z'l). We have found that the problem 
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can still be mapped onto a family of L'^ single-impurity 
one-electron effective problems defined on the very same 
lattice, Zf. Each irrep corresponds to a value of a d- 
dimensional integer vector M {Mi = ... L — 1), or more 
physically, to a value for the total momentum K = ^M. 
The only difference between irreps lies in the value of the 
hopping term which depends explicitly on M. More pre- 
cisely, for each irrep labelled by K, Eq.([l^) now reads 



[10] Identity obtained by making a series expansion of opera- 
tor {H — z)~^ in terms of V and using the fact that V is 
a projector. 



L-1 



L-1 

E 



1 



Ld ^ ^ (K) 

rii=0 nd=0 - e,7 



1 

U 



(21) 



where the non-interacting energies e^^^ are given by 



-4t cos iKj2) cos {h + Ki/2) (22) 



with K = and k — ^n. In contrast with hyper- 

cubcs, the sum in Eq.(^l|) runs over i"* values and cannot 
be further reduced. Note that, since the effective hopping 
term varies with the irrep (and can even vanish), both 
localized and resonant states (irreps by irreps) may be 
present. Finally, it is worth noticing that in the d = 1 
case the Bethe ansatz equations for two particles |^] : 



tan 



kiL 



U 



2t(sin ki — sin fc2) 



(23) 



with E = -2t(cos ki + cos ^2) and ki + k2 = (M = 

. . . L — 1) can be recovered from Eq.(^l|) after simple 
but tedious algebra. 
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